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We show that the concept of topological order, introduced to describe ordered quantum systems which cannot 
be classified by broken symmetries, also applies to classical systems. We discuss some of the fundamental prop- 
erties of this type of classical order, and propose how to expose it via a generalized topological entropy. Starting 
from a specific example, we show how to use (quantum) pure state density matrices to construct corresponding 
(classical) thermally mixed ones that retain precisely half of the original topological entropy, a result that we 
generalize to a whole class of systems. 



I. INTRODUCTION 

Over the past two decades, the notion of topological order 
has been developed to describe quantum systems exhibiting 
exotic properties such as a ground state (GS) degeneracy that 
cannot be lifted by any local perturbations'-^ and fractional- 
ized degrees of freedom.^ These systems have been shown to 
exhibit a type of order that defies the canonical classification a 
la Landau-Ginzburg, and that was named topological order. 
Levin and Wen, and Kitaev and Preskill recently proposed the 
idea that an appropriately defined topological entropy could 
be used to measure the presence of topological order in the 
GS of a quantum system.^'^ Such entropy can be obtained as a 
linear combination of Von Neumann entanglement entropies 
of different bipartitions of the system, 

'5'topo = 1™ i^^lA + S2A + S3A ~ S4A) J (1) 

r,R^ 00 

aimed at removing all bulk and surface terms to uncover the 
sole topological contribution. A particular choice of the four 
bipartitions, as used in Ref.|5l is illustrated in Fig.[Tl 
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FIG. 1 : Illustration of the four bipartitions used to compute the topo- 
logical entropy in Ref. 5. 

Despite all the efforts to understand topological order so far, 
the research focused exclusively on zero-temperature quan- 
tum systems, and the possibility of realizing topological order 
in classical systems has been left open.^ 

In this paper we explicitly show that the concept of topolog- 
ical order is not a pure zero temperature property but it applies 
also to classical systems. Starting from a known example of 
quantum topological order in a toric code, the Kitaev model, 
we show how an appropriate coupling to a thermal bath can, 
in a certain limit of the coupling constants, bring the system 
into a totally mixed state that exhibits a non-vanishing topo- 
logical entropy, therefore suggesting that topological order is 
capable of surviving thermal dephasing. We discuss one of 
the fundamental properties of classical topologically ordered 



systems, namely the peculiar structure of phase space divided 
into sectors that are connected exclusively by extensive re- 
arrangements of the microscopic degrees of freedom. These 
sectors are not distinguishable under any local measurement, 
and therefore lead to what ought to be called topological er- 
godicity breaking. A display of glassy behavior is expected 
to accompany classical topological order (and possibly vice 
versa). We briefly comment on other properties of classical 
topological order, although a thorough characterization is be- 
yond the scope of the present paper. Finally, we reformulate 
the definition of topological entropy by using, instead of the 
Von Neumann entropy, a boundary entropy that is based on 
the mutual information entropy defined in information theory. 
This allows for a natural and more symmetric extension to 
classical as well as quantum systems. We then generalize the 
results on the Kitaev model to a whole class of quantum sys- 
tems with GS wavefunctions whose amplitudes are positive 
and factorizable when expressed in terms of some (local) mi- 
croscopic degrees of freedom. We argue that a loss of half 
the topological entropy between pure states and correspond- 
ing mixed states occurs whenever such a correspondence can 
be established. 



II. TOPOLOGICAL ENTROPY IN THE KITAEV MODEL 

Consider the Kitaev model in Ref. 8. It consists of spin- 
1/2 degrees of freedom on the bonds of a square lattice with 
periodic boundary conditions. The Hamiltonian of the model 
can be written in terms of plaquette (p) and star (s) operators 
as 

H = -A^^^-As^Sp (2) 

s p 

where A^, A^ are real, positive parameters; = Yijep ^1 
where j labels the four spins belonging to plaquette p, and 
= Y[j£s '^^ where j labels the four spins belonging to the 
star centered on vertex s. 

Using the fact that all and Bp operators have eigenval- 
ues ±1, that they all commute, and that they are subject to 
the constraint Y[s = 1 = Yip ^p' oi^^ show that the 
GS of the model is degenerate and the dimension of the GS 
manifold is precisely 2^. Furthermore, one can prove that this 
degeneracy is not lifted by any local perturbation in the ther- 
modynamic limit, and the system exhibits topological order. 
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Let us first re-derive the topological entropy of the model 
with a new approach that can be immediately applied to both 
classical and quantum systems. We are going to work within 
one specific topological sector This choice will not affect 
the results as Hamma et al. showed that all sectors exhibit 
the same entanglement entropy.— In the z-basis for the spins, 
where all the operators are diagonal, the GS is given by the 
equal amplitude superposition of all states where the product 
of the four components around each plaquette is positive. 
Notice that the action of a star operator in the z-basis does 
not violate the plaquette constraint above. One can show that 
any state in the equal amplitude superposition (within a single 
topological sector) is uniquely specified by a product g of star 
operations required to obtain it from a given reference state 
|0) that satisfies as well the plaquette constraint (and lies in 
the same topological sector). The GS can thus be written as 



I*) 



1 



|G|i/; 



(3) 



see 



where G is the (Abelian) group of all possible products of star 
operators, and \G\ = 2"/^^^ is the number of elements in the 
group (n being the total number of spin degrees of freedom). 
At zero temperature, the system is described by the density 
matrix 



(4) 



As elegantly shown in Ref. |9| by Hamma et al, the Von 
Neumann entropy of any bipartition (A, B) of the system pre- 
pared in such GS is given by 



Sa = -Tr(pAlog2 Pa) 

= log2|G'| -log2(d^ds) 



(5) 



Here is obtained by taking the partial trace over subsystem 
B (p^ = Trgp) and (d^) is the number of elements in 
the subgroup C G (Gg C G) of transformations acting 
solely on A (B) and leaving B {A) unchanged. Moreover, 
Hamma and collaborators showed that whenever subsystem 
A is connected, then — 2^-* and dg — 2^^, where 
and are the number of single star operators acting solely 
on A and B respectively. 

Since we are interested in computing the topological en- 
tropy ([T]), we need to extend the above results to the case of 
bipartitions involving multiple disjoint components. Let us 
consider first the case when the set A is composed by two dis- 
joint regions Ai and A2, each of them connected, as in biparti- 
tion 4 in Fig.[T] While all the operators in G that act solely on 
A can be obtained as products of star operators acting solely 
on A, the converse is no longer true for subsystem B. In fact, 
the product of all star operators acting on Ai (i.e., both those 
stars acting solely on spins in Ai and those boundary stars act- 
ing on spins in A^ and on spins in B simultaneously) flips all 

in A-^ two times, therefore leaving them unchanged, while 
it has a non-trivial action on B. Namely, this product gives 
rise to an operator that flips a string of spins in B along the 



boundary of A^^. Similarly if we consider the subsystem A2, 
and construct the product of all star operators acting on 
The composition of these two operators gives the product of 
all star operators acting on A (i.e., both those acting solely on 
spins in A and those at the boundary acting simultaneously 
on A and on B), which is equivalent to the product of all star 
operators acting solely on B. Therefore, out of the two ad- 
ditional operators, only one of them is a new operator on B 
that cannot be obtained as a product of star operators acting 
only on B. In this specific example then, c?^ — 2^ a while 

dg = 2^B+\ 

In the generic case of A — U . . . U and B = 

-Si U . . . U B^^, with A^... A^^ and B^ . . . B^^ disjoint 
and connected, we obtain 



d^ = 2 



dB = 2' 



(6) 



Finally, we can compute the topological entropy of the sys- 
tem prepared in its GS (O using Eq. ([T]) and the four different 
bipartitions defined in Fig.[Tl From Eqs. (|5j and (|6]l, and from 
the fact that 



^lA 



^3.A 



'■4A 



= 2 



( and all others = 1), 



we obtain S,,^ 



^topo — 2 = log2(£' ), where D = 2 is the so- 
called topological dimension of the system, in agreement with 
the result by Levin and Wen J- 



ni. CLASSICAL TOPOLOGICAL ENTROPY 

Let us consider now what happens when the system in ques- 
tion is coupled to a thermal bath that allows for dephasing and 
thermalization within the lowest lying eigenstates of the pla- 
quette operators B^. This is the case in the limit of — + 00, 
i.e., acting as a local hard constraint, with A^/T 0. Later 
on we discuss how the physics of this hard-constrained system 
is relevant to the soft case where the coupling constants are fi- 
nite, and Ag 3> T ^ A^. In the hard-constrained regime the 
system is described by a totally mixed density matrix:— 



^E5|0)(0|5- 

' ' gee 



(7) 



Following the same arguments as for the pure state (01 in 
Ref. '9', we show that 



Pa 



\G\ 



E 9a\0a){0a\9a, 



(8) 



where G/Gg is the quotient group, and and |0^) are given 
by the generic tensor product decompositions of |0) = |0^) ® 
|0g) and g = Qa® Qb- We can then show that 



Pa = 



\G\ 



Pa^ 



(9) 
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and use the identity 

- Tr log2 Pa) = -^ lim ^ (xrp^i) , (10) 
to finally obtain the Von Neumann entropy 



log2 |G| -log2 

log2 dA 



^(pure state) 



B- 



(11) 

(12) 



We will comment more on this result below. From Eq. (fT2] i 
we can compute the topological entropy ([T]i: 



r*(mixed state) 
'-'topo 



r<(pure state) 
^topo 



rM, 

rt(pure state) _ , 
*~^topo ^ 



lim log. 



1. 



(13) 



Notice that, while the bulk and boundary contributions present 
in the classical Von Neumann entropy in Eq. ( fT2] l precisely 
cancel in the expression for the topological entropy, the topo- 
logical contributions give a non-vanishing result, therefore 
suggesting that topological order may survive thermal mixing 
(at least in 2D hard constrained systems). 

Notice also that the topological entropy of the classical 
model obtained upon suppressing the off-diagonal elements 
in the density matrix of the Kitaev model differs from the 
quantum case by a factor of 1/2, namely logj D instead of 
log2(-D^). The physical meaning of this factor is that each of 
the underlying gauge structures of the Kitaev model (the so- 
called electric and magnetic loops) is independently responsi- 
ble for half of the zero-t emp erature topological entropy, as we 
recently showed in Ref.[l2l 

What are the physical properties of such a classical topo- 
logically ordered system? First, it follows from the density 
matrix (|7J that the expectation value for any product of spins 
(z-component, as these are thought as classical now) is zero 
unless the spins belong to a closed, non-winding loop. The 
system is thus a featureless "spin liquid" in a classical way 
(much as a paramagnet). But it does have order, topological 
order, in the sense that there are different such paramagnets 
distinguished by the product of spins over a winding loop on 
the torus. When the constraint rij^p^'j = +1 is fully en- 
forced (Xg/T — > oo), classically changing the topological 
sector requires the concomitant flipping of a number of spins 
that is of the order of the linear size L of the system. The prob- 
ability for such event is of order q^, where q''="s*(f) jg jjjg 
dependent probability of the combined flipping of all the spins 
along a loop £ (equivalent to a quantum tunneling event). This 
leads to the breakdown of phase space into disconnected topo- 
logical sectors (broken ergodicity), and the associated time 
scales grow exponentially in the size of the system, a typi- 
cal signature of glassinessji^ Notice the topological nature of 
such ergodicity breaking process, in that the disconnected sec- 
tors cannot be distinguished by any local measurements. 

Other properties observed in quantum topologically or- 
dered systems, such as fractionalized excitations, lead to less 
intuitive classical counterparts to be found, for example, in 
the behavior of defects. In our classical version of the Kitaev 
model they appear in pairs connected by fluctuating strings. 



and they are likely to give rise to unusual response and re- 
laxation processes. While the topological entropy vanishes 
identically, in the thermodynamic limit, as soon as the con- 
straint is softened, '^''^ the topological nature of the defects 
will persist at small enough temperatures (A^ ^ T ^ A^^), 
very much as the fractionalized excitations do in the quantum 
case (A^, A^ ^ T). In fact, at such low temperatures the de- 
fects are sparse and the constraint is satisfied everywhere in 
between. Thus the defect creation, motion and annihilation 
is locally dictated by the same "vacuum" structure that would 
lead to topological order, were it enforced throughout the sys- 
tem. 



IV. CONSTRUCTING A "BOUNDARY ENTROPY" 

From Eq. ( fT2] l. one can see that the Von Neumann entropy 
of a classical system is not a good measure of the boundary 
entropy of a bipartition of the system, because of the explicit 
dependence on the bulk entropy of one of the two partitions. 
One can obviate the problem by using the mutual information 
entropy of a bipartition, defining the boundary entropy to be 



-^boundary 



(14) 



Clearly S'l^^j^^j^^ = Sj^ = Sg when the system is prepai'ed in 
a pure state. In the case of a thermally mixed state, all bulk 
entropy contributions cancel in S'jjpy^j^jy, contrarily to the Von 
Neumann entropies S"^ and Sg, and only boundary terms are 
retained. 

We propose here to use .Sijounciaiy 1" Eq. ( fT4] i as an alternative 
definition of the Von Neumann entanglement entropy, appli- 
cable to quantum as well as classical systems. With this defi- 
nition, Eqs. Q and ( fT2b imply that the entanglement entropy 
stored in the boundary of a bipartition of a system prepared 
in the pure equal amplitude superposition state ^ is twice as 
large as the entropy stored in the classical counterpart (|7]): 



'{mixed state) 
boundary 



1 o(pure 



state) 



boundary 



(15) 



At least half the entanglement entropy in the Kitaev model has 
actually a classical origin''*! 

Accordingly, we redefine S^^^^ in Eq. ([T]l using Eq. ( fT4l l 



-^topo 



lim 



(1) 



:.(2) 



(3) 



(4) A 



^, y boundary ' boundary ' boundary boundary J 

(16) 

For the Kitaev model we immediately obtain 



o(inixed state) ^ C' 

'-'topo "^^-'topo 



1 X 



-(pure state) 



(17) 



in agreement with our direct calculation in Eq. ( fTST l. We argue 
that Eq. ( fT6] l is the proper reformulation of the topological 
entropy by Levin and Wen in order to extend it to classical 
and quantum systems alike. Together with Eq. ( fTSl ). which 
we show hereafter to hold for a wide class of GS wavefunc- 
tions, the new definition of topological entropy implies that 
also Eq. (T7\ holds true within the same class. For this type 
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of systems, at least half of the quantum topological order has 
therefore a classical origin. (Whether part of the remaining 
half of the topological entropy could be further ascribed to 
classical correlations, in the sense of Ref. 14, is an interesting 
problem.) 



V. BEYOND THE KITAEV MODEL 



All the calculations, as well as Eqs. (I3ll5l ITIfTSIi. generalize 
straightforwardly to any quantum system whose GS is given 
by the equal amplitude superposition of a given set of states.— 

On the other hand, extending our results to generic wave- 
functions \^!) ~ J2geG ^is) 9\^) highly non- 
trivial because g = g^^ ^ does not necessarily imply that 
a{g) = a((7^) a{gg)- Supporting results in this direction 
come from considering the case of non-negative, local wave- 
functions, as defined by the condition that in the continuum 
limit aig) ^ a(</)) = 0-^(0^) a^(0^) a^(0-^), with o^(0^) 
depending only on the boundary d between subsystems A and 
B where the two field configurations (j)^ and (j)^ match (and 
equal 0^).'^ This is the case, for example, of scale invariant 
GS wavefunctions, such as those of systems at a quantum crit- 
ical point. 

Consider the GS wavefunction(al) in the continuum limit 



V Z . 



(18) 



where Z = J D(f)e~'^^^'f'K and E{(j)) satisfies E{<j)) = 
E'^{(l)^) + E^{(f)^) + E^{(l)'^,(j)^)withE^{(j)^,(t)^) depen- 
dent only on the boundary d between subsystem A and sub- 
system B, for any bipartition {A, B). 

As discussed in Ref.[l7|, one can start from the correspond- 
ing pure state density matrix 



P = 



(19) 



and obtain the elements of the reduced density matrix = 

{^t\PA\<l>2) = 



-/3[£;-'^(0^)+£;«(0«)+_E°(0^,0«)]/2 



(20) 



One can then evaluate 



1—1 

n 



(21) 



where by construction = (pi and ipn+i = The term 
nr=i ^i'f'i ^ 't'i+i ) is the boundary constraint arising from the 



product of with itself (which enforces the boundaries to 
all match, because (pf matches which then matches (f>f^i, 
which in turn matches (pf^i, and so on). If we factorize the 
integral over Y[7=i ^'^i Y[7=i D(j)fD(t)fD(j)f, we can sim- 
plify the integrals over boundaries and the product of bound- 
ary delta functions to obtain 

(22) 

where we explicitly used the assumption that all the 
E^ ) terms depend solely on the respective bound- 

aries 0f . Now that the delta functions are no longer present, 
we can address the remaining integrals as follows. The con- 
dition that the n fields 0^ and the n fields agree at the 
boundary is equivalent to having 1 replica free to take any 
value and the remaining n — 1 to be pinned, with fixed (Dirich- 
let) boundary conditions. One thus obtainsi^ 



- I D<l)^D(f)'^D<j)^ e 

z 



-PnE^-l3E'^(<j,'^)-l3E^{(t>^) 



(23) 



where Z^ — Jj^ Dcj)^ ^ (<p ) gjygjj i^y (jjg integral over 
subsystem A with Dirichlet boundary conditions, and equiva- 
lently forZ|. 

Finally, using the identity ( fTOt . we can compute the Von 
Neumann entropy 



7 A yB 




Z 




yA yB 




Z 





Sa = -log2 



= -log2 

The case of the totally mixed state 



Z J ^ln2 



P = 



D(j)e 



-liE{4>) 



(24) 



(25) 



involves similar calculations, where the term 
nLi'^('/'f''^f+i) in Eq. (EB gets replaced by 

nr=i^('/'f,d)' 



<^i+i)- a result 



B \ 



Z 



z 



and 



Sa = -log^ 



yB 

z ■ Z 



(26) 



PE{4>) 



(27) 
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Notice that, while for an equal amplitude superposition the 
totally mixed state could be reached via an appropriate high- 
temperature limit (i.e., via coupling to a thermal bath), this 
is no longer true for Eqs. ( fT9] l and (IZST i, which should be 
interpreted only as a recipe to construct an associated classical 
system, given the quantum one. 

In conclusion, we obtain 

^(pure state) ^ pp^ ^ pp^ ^ pp^^^ ^ ^^J^a^ 
^(^ixedstate) ^ pp^ _ pp^^^ ^ ^^^pB^ ^ ^^^pA^ ^ 

Where l3Fj, - - log2 Z^, PP^ = - loga Z%, and /^F^^b - 
-log2 Z. 

For the boundary entropy defined in Eq. <fT4l > we find that 



oCmixed state) 
boundary 

__ 1 

" 2 





+ A 

In 2 




ln2 




+ A 

ln2 




ln2 


In 2 


l^pAuE 







AUB 



- (^/3F^+;3Fb-/3F^ub + i^(^^ 

Q(pure state) 
2 boundary 



(28) 



Thus, once again, for the boundary entropy defined in Eq. (fT4l i 
we find that Eq. dTSl ) (and therefore Eq. ([TT] )) 5f//Z holds. For 
any quantum system whose GS wavefunction is local and 
with positive amplitudes, there exists an associated classical 
system which exhibits precisely half the topological entropy. 
The configurations of the classical system have as Boltzmann 
weights the squares of the amplitudes in the associated quan- 
tum GS. 



the definition of topological entropy^'^ that applies to classical 
and quantum systems alike. We use this topological entropy 
to identify topologically ordered phases at the classical level, 
and we discuss some of the properties of such classical sys- 
tems. We also show how to construct thermally mixed (classi- 
cal) density matrices that retain precisely half the topological 
entropy of associated pure state density matrices for a special 
class of quantum systems. Our results imply that for quan- 
tum systems whose GS wavefunction in a given basis can be 
gauged so as to have non-negative (local) amplitudes, quan- 
tum topological order is present if and only if the correspond- 
ing classical system obtained by removing all off-diagonal 
matrix elements of the quantum density matrix is also topo- 
logically ordered as defined above. 

The result that topological order applies to classical sys- 
tems could possibly provide a way to classify classical orders 
without an obvious order parameter We end with a specu- 
lative note on the possibility that this might be the case in 
glassy systems. No order parameter can be constructed in 
glasses from equal-time correlation functions (as opposed to 
the two-time Edwards-Anderson order parameter) of physical 
degrees of freedom (as opposed to replicated variables). This 
fact is suggestive that, much as in quantum topologically or- 
dered systems like spin liquids, no local order parameter can 
detect the underlying order of the glassy states, and the hid- 
den order could indeed be topological. In this case, our finding 
that topological order and topological entropy can be defined 
in classical systems may have implications to understanding 
the physics of glassy systems. 
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